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Abstract 



Levinson's theorem for Dirac particles constraints the sum of the phase shifts 
at threshold by the total number of bound states of the Dirac equation. Re- 
cently, a stronger version of Levinson's theorem has been proven in which 
the value of the positive- and negative-energy phase shifts are separately con- 
strained by the number of bound states of an appropriate set of Schrodinger- 
like equations. In this work we elaborate on these ideas and show that the 
stronger form of Levinson's theorem relates the individual phase shifts di- 
rectly to the number of bound states of the Dirac equation having an even or 
odd number of nodes. We use a mean-field approximation to Walecka's scalar- 
vector model to illustrate this stronger form of Levinson's theorem. We show 
that the assignment of bound states to a particular phase shift should be 
done, not on the basis of the sign of the bound-state energy, but rather, in 
terms of the nodal structure (even/odd number of nodes) of the bound state. 

PACS: 11.80.-m, 11.10.Q 
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I. INTRODUCTION 



One of the most beautiful results in scattering theory was proven by Levinson in 1949 [|TJ . 
Levinson's theorem relates the number of bound states in a specific angular-momentum 
channel (Ni) to the value of the phase shift (5i) at threshold 

$,(0)=JV,7T. (1) 

In arriving at this form of Levinson's theorem the modulo- 7r ambiguity in the definition of 
the phase shift has been resolved by demanding that the phase shift approaches zero (and 
not only a multiple of n) in the high-energy limit. For simplicity, we have assumed, and will 
assume hereafter, the absence of zero-energy resonances. 

Finding a generalization of Levinson's theorem to Dirac particles, however, has proven 
to be a subtle task. Initially, it was believed that the generalization of Levinson's theorem 
related the difference of positive- and negative-energy bound states to the corresponding 
difference of positive- and negative-energy phase shifts at threshold @. It was also suggested 
that the modulo- tc ambiguity in the definition of the phase shift could be resolved, as in the 
nonrelativistic case, by allowing the phase shift to vanish in the high-energy (E —>■ ±00) 
limit. These results, however, were later found to be incorrect and it was only until 1985 
that a correct generalization of Levinson's theorem to Dirac particles was found by Ma and 
Ni ||. They showed that the correct form of Levinson's theorem for Dirac particles relates 
the sum (and not the difference) of the phase shifts at threshold to the total number of 
bound states [Q] 

5 K (E = +M)+5 K (E = -M) =N K n. (2) 

In the above equation k is the generalized angular momentum channel (j — \k\ — 1/2) and 
N K is the total number of (positive- plus negative-energy) bound states in that channel. 
In this form of Levinson's theorem one implicitly assumes that the modulo- n ambiguity in 
the definition of the phase shift has already been resolved. However, in contrast to the 
nonrelativistic case, where one is free to define the high-energy limit of the phase shift to be 
zero the asymptotic value for the phase shift is given by |||7|]. 

lim 6 K (E) = T / V(r)dr } (3) 

E^±oo Jo 

where V(r) is the timelike component of a Lorentz-vector potential. In particular, this last 
relation implies 

5 K (E ^ +00) + 5 K (E^ -00) = 0, (4) 

and indicates that it is only the sum, but not the individual phase shifts, that vanishes in 
the high-energy limit. Yet, Eq. @ establishes that the asymptotic behavior of the of the 
individual (positive- and negative-energy) phase shifts is known. In contrast, the behavior 
of the individual phase shifts at threshold is (as opposed to their sum) not constrained by 
Levinson's theorem. Recently, however, Poliatzky was able to prove a stronger version of 
Levinson's theorem by showing that the individual phase shifts at threshold are related 
to the number of bound states of an appropriate set of Schrodinger-like equations which 
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coincide with the Dirac equation in the limit of zero momentum ||. In contrast to previous 
expectations 0, however, the phase shift at threshold with a particular sign of the energy 
does not coincide (in general) with the number of bound states of the Dirac equation having 
the same sign of the energy, i.e., 

8 K (E = ±M) ^ iVf )tt . (5) 

As we will show later, this "naive" expectation is upset whenever the potential becomes 
sufficiently strong to cause the binding energy of certain states to exceed the rest mass of 
the particle. 

In the present work we show, using some of the ideas developed by Poliatzky |§, that 
the behavior of the phase shift at threshold with a particular sign of the energy is directly 
related to the number of bound states of the Dirac equation having a specific nodal structure. 
We will establish that the assignment of bound states to a particular phase shift should be 
done, not on the basis of the sign of the bound-state energy but, instead, by determining 
if the combined number of nodes in the upper and lower components of the bound-state 
wavefunction is even or odd. As we shall see, this criterion is inextricably linked to the 
orthogonality and boundary conditions satisfied by the eigenstates of the Dirac equation. 

Our paper has been organized as follows. In Sec. [TT] we review most of the formalism 
needed to prove the stronger form of Levinson's theorem. Particular emphasis is placed on 
the boundary conditions satisfied by eigenstates of the Dirac equation which are ultimately 
responsible for the nodal structure of the bound state. In Sec. |T| we illustrate Levinson's 
theorem using a mean-field approximation to the Walecka model as an example. In par- 
ticular, we show that the "naive" relation between the phase shift at threshold and the 
corresponding number of bound states with the same sign of the energy is upset when- 
ever the potential becomes strong enough to generate binding energies that exceed the rest 
mass of the particle. Finally, in Sec. [TV] we offer our conclusions and speculate on future 
applications of Levinson's theorem. 



II. FORMALISM 



In the presence of a spherically symmetric potential the eigenstates of the Dirac equation 
can be classified according to a generalized angular momentum k and can be written in a 
two component representation 



9ek\T) 



r 



(6) 



with the upper and lower components expressed in terms of spin-spherical harmonics 

I . f k if k > 0; 



y Km (r) = (f|Z|jm) , j — \k\ — 



I 



■(1 + k) if k < 0. 



(7) 



For a particle moving in the presence of Lorentz scalar (S) and time- like vector (V) potentials 
the Dirac equation 
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a-P + f3M + pS(r) + V(r) 
reduces to a set of first order, coupled differential equations 

d k 



-T + -) 9EK{r 
dr r 



d 

dr r 



M*(r) + E*(r) 
f EK (r)= M*(r)-E*(r) 



fEn{r) 



9En{r) 



written in terms of a local mass and a local energy defined by 

M*(r) = M + Sir) , E*(r) = E — V(r) . 



(8) 

(9a) 
(9b) 

(10) 



It is a well known fact that the above set of coupled equations can be reduced to a single 
second-order differential equation identical in form to a Schrodinger equation. For example, 
one can use Eq. ( |9aT ) to express the lower component of the Dirac wavefunction in terms of 
the upper component 



fEn{r) 



1 



M*(r) +E*(r) 
9En(r) = [£(r)} 1/2 u EK (r) = 



d k 
Tr + ~r 

M*(r) + E*(r 



, + - 9ek\t) 
dr r 



1/2 



UEnir) 



(11a) 
(lib) 



which upon substitution into Eq. 
(transformed) upper component 



M + E 

yields the following Schrodinger-like equation for the 



<P_ 2 

dr 2 



k(k + 11 



- U eS (E,K;r) 



u EK (r) = ; (f = E 2 - M ; 



(12) 



1 /2 

The transformation factor [£(r)] relating the Dirac upper component (g EK ) to the 
Schrodinger-like component (u Ek ) was introduced in order to remove all terms linear in 
the derivative of g En . The effective potential, which has now acquired an energy dependence 
due to the transformation, 



U eS {E, k; r) = U C {E; r) - kU so (E; r) + U D (E; r) 



(13) 



is defined in terms of central, spin-orbit, and Darwin contributions that are, respectively, 
given by 



U c {E-r) = 2M 
U so (E;r) = -- 



S(r) + ^V(r) ) + 



'S 2 (r)-V 2 {rY 
2M 



1 dtjjr) 
£(r) dr 



U D (E;r) 



1 d£(r)~ 


2 1 


[ 1 d 2 £(r)l 


^(r) dr 


2 





(14a) 
(14b) 

(14c) 
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In particular, the energy dependence displayed by the (central) potential now makes evi- 
dent the high-energy limit of the phase shifts [Eq. (HD) through the well known asymptotic 
behavior of the Schrodinger phase shift [EM 



lim 5 K (E) = / V cS (r)dr^ T \ 

^±00 p Jo Jo 



V(r) dr . 



(15) 



Thus, the asymptotic behavior of the phase shift can ultimately be traced to the Lorentz 
transformation properties of the vector potential. In addition, this relation indicates that 
the high-energy limit of the phase shift is unaffected by the presence of a Lorentz-scalar 
potential. Moreover, this finding contradicts the nonrelativistic notion that a distortion-free 
description of the projectile should be appropriate at large enough energies. 

In reducing the original couple set of equations to an effective Schrodinger-like equation 
one could have, alternatively, chosen to eliminate the upper component of the Dirac wave- 
function in favor of the lower component. This can be achieved, because of the structure of 
the equations [see Eq. @], by simply changing the sign of k and of all timelike quantities 
(i.e., E — > —E and V — > —V). These modifications give rise to the following relations 



9ek{t) 



M*(r) - E*(r) 
1/2 



fEnir) = £(r) v E n{r) 



~M*(r)-E*(r)~ 
M-E 



1/2 



VEn{r) 



(16a) 
(16b) 



and generate the corresponding Schrodinger-like equation for the (transformed) lower com- 
ponent 



dr 2 



R(R + 1) 



- U eS (E, k; r) 



VEkVT) 







(17) 



In the above equation we have defined R = — k and introduced the effective potential U e R 
that is obtained from U e s by changing the sign of all timelike quantities. 

The system of equations Eqs. flTT|) and (|T^), or equivalently Eqs. fllBD and ([1?]), with 
the appropriate boundary conditions, are equivalent to the original set of coupled equations 
[Eq. (|9])]. Unfortunately, the nonrelativistic form of Levinson's theorem can not be applied 
to these equations. The difficulty stems from the fact that the theorem is not valid once the 
effective potential has acquired an energy dependence. Fortunately, Poliatzky has recently 
been able to show how to overcome this difficulty [§]. He has proven that for an equivalent 
set of Schrodinger-like equations, 



d 2 9 k(k+1) tt /rl 
d^+P I* " U *( E = +M ' K ' r ) 

d 2 



dr 2 ^ ^ 



r(r + 1) _ 



U eS (E=-M, R;r) 



uek{t) = , 

VEn{r) = , 



(18a) 
(18b) 



obtained from replacing the energy- dependence of the effective potential by its value at 
threshold, the nonrelativistic form of Levinson's theorem can indeed be applied ||. This 
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stronger form of Levinson's theorem constrains, not only the sum of the phase shifts [Eq. (Q, 
but, in addition, the behavior of the individual phase shifts at threshold 

5 K (E = +M) = 5 { l +) (p = 0) = n{ +) vr , (19a) 
8 K (E=-M) = 5\~\p = 0) = nf ] Ti . (19b) 

In the above equations the orbital angular momentum I and I are related, respectively, to 
k and k through Eq. ([?[). In addition, n| and n[ ^ are the number of (nonrelativistic) 
bound states supported by the energy-independent potentials {U e s and U e $ respectively). 
The equality between phase shifts is a direct consequence of the equivalence of the Dirac 
equation to Eq. ( |18a| ) and Eq. (|18b| ) at the corresponding thresholds. The relation to the 
number of bound states follows directly from applying the nonrelativistic form of Levinson's 
theorem which is now valid because the potentials have become independent of energy. 

We now show that the consistency of the above set of equations enables one to relate the 
value of the phase shifts at threshold to the number of bound states of the Dirac equation 
having, either, an even or odd number of nodes. Essential to the proof is the understanding 
of the asymptotic behavior (both as r — > and as r — > oo) of the ratio of lower to upper 
components of the Dirac wavefunction. 

For small values of r the upper and lower components of the Dirac wavefunction can be, 
respectively, written as 

g EK (r) ~ Aji (p*r) ; f EK (r) ~ Bjj (p*r) , (20) 

where ji(z) = zji(z) is the Ricatti-Bessel function (i.e., spherical Bessel function times its 
argument) ||. The coefficients A, B and the wavenumber p$ are determined (after using the 
Dirac equation), by solving the secular equation 

Pq sgnft — Mq — Eq 
-Mo + Eq -pi sgnft 

which results in 

p* = V^ 2 - , (22) 



(21) 



Po 



AfJ + E, 



or 



o , 



-sgnK 



Po 



M$ - E, 



(23) 



o , 



where we have defined sgn/t = k/\k\, Mq = M*(r = 0), and Eq = E*(r — 0). 

The asymptotic (r — > oo) behavior of the wavefunction can be obtain in a similar fashion. 
In this case, however, one must distinguish between scattering states (E > +M or E < — M) 
and bound states (— M < E < +M). For scattering states the asymptotic behavior is given 
by 



9Eti{r) ~ Csin 



In 

pr - — + S K (E) 



f EK {r) ^ ^Dsin 



pr 



y + S K (E) 



(24) 
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and yields the following relations, 



p = VE 2 - M 2 , (25) 



= sgn/t ( — ^— ) , or (-) = sgnK ( — ^— ) . (20) 



.C7 ° \E + MJ ' V-D/ \E-M, 

For bound states, on the other hand, the asymptotic behavior is given by 

^(r) ~ Ce~ pr ; /^(r) ~ De~ pr , (27) 

which in turn implies 



p = VM 2 - E 2 , (28) 



°' ' " ' o, ^U-^). (29) 



C7 \M + EJ ' V-D/ \M-E 

In particular, this last relation establishes (because — M < E < +M) that the ratio of the 
two bound-state components is always negative irrespective of the value of k. 

Having established the boundary conditions satisfied by the two components of the Dirac 
equation we now turn to the set of equations that have determined the threshold behavior 
of the phase shifts. We start by inspecting Eq. (|llb| ) which defined the relation between 
qek and ue k and ultimately determined the value of S K (E = +M). The fact that both gE K 
and uek satisfy real differential equations with real boundary conditions ||, demands that 
the transformation factor £(r) = [M*(r) + E*(r)]/[M + E] be positive for all values of r 
and in particular positive for r = 0. This implies that bound states of the Dirac equation 
satisfying this condition, and therefore being "counted" by 8 K (E= +M), have the sign of 
the ratio of the two Dirac components at the origin determined solely by the sign of k, i.e., 
from Eq. (p3|) we obtain 



sgn 



sgn I — = sgn/t 



r— >0 



Since, in addition, the ratio of the two Dirac components is always negative at large separa- 
tions [see Eq. (|29|) 1, these relations constrain the combined number of nodes in the bound- 
state wavefunction (i.e., the number of nodes in qek plus the number of nodes in to be 
even for k<0 and odd for k>0. In the case of k>0, for example, the ratio is positive near 
the origin and the lower component must change sign an odd number of times more than 
the upper component in order for the ratio to become negative at large separations. 

Relating the behavior of S K (E = —M) to the number of bound states with the compli- 
mentary nodal structure now follows from Eq. (|16b| ). Again, the fact that both f EK and Ve k 
satisfy real differential equations with real boundary conditions forces the transformation 
factor £(r) = [M*(r) — E*(r)]/[M — E] to be positive for all values of r. Consequently, the 
behavior of the ratio of Dirac wavefunctions at the origin is now determined, not by the sign 
of k but, instead, by the sign of R = —k, i.e., 
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sgn 



9eJ/) 



= sgn (— \ = -sg: 



Consequently, bound states of the Dirac equation satisfying the above relation, and hence 
being counted by S K (E = — M), must have a combined number of nodes that is now even for 
k>0 and odd for k<0. Thus, the boundary conditions have proved essential in establishing 
the connection between the threshold behavior of the phase shifts and the number of bound 
states of the Dirac equation with a specific nodal structure. 

We can now formulate the stronger version of Levinson's theorem for a Dirac particle 
as follows. If the modulo- 7r ambiguity in the definition of the phase shifts is resolved by 
demanding the fulfillment of the high-energy relations 



and if N^(Njf') represents the number of bound states of the Dirac equation having an 
even(odd) number of nodes, then the phase shifts at threshold are given by 



So far we have ignored the possible existence of zero-energy resonances. The presence of 
zero-energy resonances can only occur in channels having, either, k = — 1 or k = +1 and 
their existence modify the above relations by the addition of a factor of tt/2 in the value of 
S K (E = +M) or S K (E = -M) respectively . 

Deeply rooted in the proof of Levinson's theorem is the completeness of eigenstates of 
a Hermitian Hamiltonian. The completeness property establishes that the appearance of 
a bound state must always be accompanied by the corresponding disappearance of a state 
in the continuum. Alternatively, the existence of an attractive potential does not alter the 
(infinite) number of states but simply "pulls down" some scattering states into the bound- 
state region @^J]. The phase shift at threshold may, thus, be regarded as a bound-state 
"counter" that triggers every time a scattering state is pulled down into the bound-state 
region. 

Examined against this background, and with the added constraint imposed by orthog- 
onality, it may not be surprising that the phase shifts at the two different thresholds 
(E = ±M) are related to the bound states of the Dirac equation with a complimentary 
nodal structure. As the potential becomes strong enough to support the existence of bound 
states, the phase shifts at threshold will "trigger" every time a scattering state is pull down 
into the bound-state region. States moving into the bound-state region from different ends 
of the continuum, however, must necessarily have a different nodal structure. This is a 
consequence of the orthogonality of states with different energies. For example, consider the 
first state (i.e., the most deeply bound) that has been pulled down from the negative-energy 
continuum and the corresponding one from the positive-energy continuum. These are the 
bound states that are expected to have the fewest number of nodes [see, e.g., Fig. [|. These 
two states, however, can not be simultaneously free of nodes; otherwise they could not be 
orthogonal and still satisfy the asymptotic condition [Eq. (p9[) ]. Therefore, all states that 
have been pulled down from the negative-energy continuum must have a complimentary 




(30) 




(31) 



S 



nodal structure to those states that have originated in the positive-energy continuum. Fur- 
thermore, this nodal structure is preserved even when the potential is strong enough to cause 
the binding energy of some of the states to exceed the rest mass of the particle. Therefore, 
it is the nodal structure (even or odd number of nodes) and not the energy (positive or 
negative) that determines the behavior of the phase shift at threshold. 

Implicit in the above derivation of the stronger form of Levinson's theorem is the assump- 
tion that the wavenumber in the interior (j)q) is real. This fact conforms to the intuitive 
notion that a bound-state wavefunction oscillates in the interior and falls off exponentially 
in the exterior. In particular, this implies that (Mq+Eq) and (Mq — Eq) must have opposite 
signs; a fact that follows directly from Eq. (p2|), i.e., 

- pf = (M* 2 - E* 2 ) = (M* + El) (M g * - E*) < . (32) 

Consequently, all bound states can be classified according to, either, (Mq + Eq) > [those 
related to 5 K (+M)] or according to (Mq - Eq) > [those related to 5 K (-M)]. 

There are, however, some very interesting bound states for which the interior wavenum- 
ber, as defined by Eq. (0), is purely imaginary. These bound states, which appear in the 
context of strong-coupling field theories, display exponential behavior, not only at large dis- 
tances, but, in addition, in the (interior) region where one customarily observes oscillating 
behavior ||,[10) . These shell-like states are driven by a strong scalar potential that generates 
a strong binding energy and, thus, a rapid variation of the wavefunction. The existence of 
nodes (which leads to an increase in the kinetic energy), of an angular momentum barrier, or 
of an additional vector potential are all factors that may destroy the shell-like nature of these 
states @. These node-free states are, thus, believed to exist only in channels with k = ±1 
and, as we show later, do not seem to violate the stronger form of Levinson's theorem. 

We now proceed to illustrate some of the ideas developed so far by means of some 
examples. 



III. RESULTS 

To illustrate Levinson's theorem we use a mean-field approximation to the Walecka 



model PT] , im . The Walecka model is a relativistic quantum field theory of nucleons in- 
teracting via the exchange of scalar and vector mesons. The ground-state of the system 
is obtained from solving self-consistently a set of mean-field (Hartree) equations The 



outcome of the calculation is a set of single-particle Dirac orbitals together with scalar and 
(timelike) vector mean fields. In Fig. [I] we show the self-consistently determined mean fields 
for 40 Ca. The potentials are strong (~ M/2) and opposite in sign in order to reproduce the 
weak-binding energy but strong spin-orbit splitting characteristic of single-particle nucleon 
(i.e., positive-energy) states. 

We have used these spherically-symmetric potentials to generate all bound states of the 
Dirac equation having k = — 1. In Fig. [2] we display the single-particle spectrum generated 
from these mean-field potentials by scaling, both scalar and vector fields, by factors of 
0.6, 1.0 (unsealed) and 1.4 respectively. The two numbers enclosed in brackets correspond, 
respectively, to the number of nodes in the upper (qek) an d lower (/ek) components of the 
bound-state wavefunction. The fact that the vector field, but not the scalar, changes sign 
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under charge conjugation leads to spectroscopy of negative-energy states that is driven, in 
contrast to the positive-energy states, by a very strong central attraction and weak spin-orbit 
splitting. In particular, for A = 1.4, the state with the least number of odd nodes (and, thus, 
"originally" in the negative-energy continuum) has crossed into the positive-energy region. 
Figures |3], f|, and [5], display the energy dependence of the positive- and negative-energy 
phase shifts and illustrate the stronger form of Levinson's theorem (see also Table |). These 
figures show how the phase shifts at threshold act as bound-state counters once the modulo- 
7i ambiguity in the definition of the phase shift has been resolved according to Eq. (|3|). For 
the first two cases (A=0.6, 1.0) the potential (although quite strong) is not strong enough to 
bind any Dirac orbital by more than the rest mass of the nucleon. Hence, the classification 
of bound states can be done according to their nodal structure or, equivalently, according 
to the sign of the energy, i.e., 

5 K (E = +M) = N^ = iV(+) , (33a) 
6 K (E = -M) = N^ o) = JVH . (33b) 

For A = 1.4, however, the binding energy of one bound-state orbital (the one with one single 
node) exceeds the rest mass of the nucleon and it becomes inappropriate to count bound 
states according to the sign of the energy (see Table |), i.e., 

5 K (E = -M)=Nj°^NiK (34) 

So far we have only studied bound-state orbitals that display (the traditional) oscillating 
behavior inside the region of the potential and exponential falloff in the exterior. We now 
turn to the very interesting case of shell-like bound states. For simplicity, we consider a 
relativistic Lorentz-scalar square well. The width (c) of the potential is fixed at cM = 20. 
The strength of the potential, on the other hand, is varied from < |So|/M< 2. Fig. |6| shows 
the K——1 single-particle spectrum along with the threshold behavior of the positive-energy 
phase shift as a function of the scalar strength. Bound states may be obtained (graphically) 
by finding the intersection of a circle of radius R 2 = (M 2 — M* 2 )c 2 with tangent-like lines. 
In particular, this expression reduces, in the weak-coupling limit (\So\/M « 1), to the 
well-know nonrelativistic relation 



R 2 



M 2 - M* 2 



c 2 ^2M\S \c 2 . (35) 



Hence the behavior of the spectrum is, at least for \Sq\/M < 1 easily understood. As 
the strength of the potential increases the radius R increases accordingly and gives rise to 
the formation of bound states. All these states are characterized by oscillating behavior 
inside the region of the well and the appearance of any bound state is accompanied by a 
corresponding increase in the value of the phase shift at threshold. As \Sq\/M = 1 (and 
M* = 0) the radius attains its maximum value and leads to the maximum number of bound 
states that can be supported by this (cM = 20) potential. As the strength of the potential 
increases even further, however, the radius R decreases leading to a weaker binding and 
ultimately to the disappearance of all but one bound state. This interesting dynamics is 
nicely reflected in the threshold behavior of the positive-energy phase shift. In fact, since for 
pure scalar potentials the binding energy never exceeds the rest mass of the particle, even 
the naive expectation for Levinson's theorem [Eq. (0)] is satisfied. 
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The sole remaining bound state is characterized by having a negative effective mass 
inside the region of the potential and satisfies < E < |M*| = — M*. In particular, the 
effective wavenumber p$, as defined in Eq. (|2^), is purely imaginary and leads, as shown 
in Fig. [5], to the formation of the shell-like state. Since the nodal structure of this state 
remains unchanged during the oscillating-like to exponential-like transition, the existence of 
this shell-like state is still properly accounted for by the value of the phase shift at threshold. 

Notice that for this shell-like state the transformation factor £(r), relating the upper 
component of the Dirac wavefunction to the corresponding Schrodinger-like wavefunction 
[Eq. ( p. lb )], is purely imaginary in the interior and purely real in the exterior. This suggests 
that, although it should be possible to obtain this shell-like solution from solving the equiv- 
alent Schrodinger-like equation, one should exercise care in enforcing these very peculiar set 
of boundary conditions. 



IV. CONCLUSIONS 

Levinson's theorem for Dirac particles relates the sum of the positive- and negative- 
energy phase shifts at threshold to the total number of bound states of the Dirac equation. 
Recently, Poliatzky was able to prove a stronger version of Levinson's theorem in which the 
value of each individual phase shift at threshold is related to the number of bound states of a 
pair of Schrodinger-like equations which coincide with the Dirac equation at zero momentum. 
In this paper we have elaborated on some of these ideas and have shown that the stronger 
version of Levinson's theorem relates the value of the Dirac phase shift at threshold to the 
number of bound states of the Dirac equation with a specific nodal structure (i.e., even or 
odd number of nodes). 

We have shown that the simple picture of Levinson's theorem developed in the nonrela- 
tivistic context and based on the completeness relation is preserved in the relativistic case. 
This picture suggests that the existence of an attractive potential does not alter the number 
of states but simply pulls down some scattering states into the bound-state region. The 
merit of Levinson's theorem is to identify the phase shift at threshold as the bound-state 
counter. 

We have illustrated the stronger form of Levinson's theorem using a mean-field approx- 
imation to the Walecka model. We have explicitly shown that the "naive" nonrelativistic 
generalization of Levinson's theorem, namely, 5 K (E = ±M) = N^ir, is violated as soon as 
the binding energy exceeds the rest mass of the particle. The nodal structure of the bound 
states, however, is robust and remains unchanged even when the qualitative structure of the 
bound state is modified (e.g., shell-like solutions). 

One essential feature that must be addressed in the study of Levinson's theorem for 
Dirac particles is the asymptotic behavior of the phase shift [Eq. (|3])]. The fact that the 
phase shift remains finite even at very large energies is a relativistic effect that might lead to 
interesting consequences. For example, many coincidence (e, e'p) experiments have already 
been proposed at CEBAF in the hope of identifying novel behavior in the propagation 
of nucleons through the nuclear medium (color transparency). To date, most relativistic 
analyses of the (e, e'p) reaction employ an ejectile wavefunction obtained from solving a 
Dirac equation having Lorentz scalar and (timelike) vector potentials. Because most of 
the analyses conducted so far have been restricted to low-energy ejectiles, the consequences 



11 



of the high-energy behavior of the phase shift have not yet been fully explored. As the 
energy of the ejectile increases, as seems to be required for the onset of color transparency, 
it will become essential to understand the relativistic effects associated with the high-energy 
propagation of the ejectile. 
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FIGURES 



FIG. 1. Self-consistent scalar and vector potentials for °Ca obtained from a mean-field ap- 
proximation to the Walecka model. 

FIG. 2. Single-particle spectrum for the k = — 1 channel in 40 Ca for different values of the 
scaling parameter A. 

FIG. 3. Positive- and negative-energy (k = — 1) phase shifts as a function of energy for A = 0.6. 
The dashes represent the asymptotic values of the phase shift (i.e., ± the integral of the vector 
potential). 

FIG. 4. Positive- and negative-energy (k = — 1) phase shifts as a function of energy for A = 1.0. 
The dashes represent the asymptotic values of the phase shift (i.e., ± the integral of the vector 
potential). 

FIG. 5. Positive- and negative-energy (k = — 1) phase shifts as a function of energy for A = 1.4. 
The dashes represent the asymptotic values of the phase shift (i.e., ± the integral of the vector 
potential). 

FIG. 6. Single-particle spectrum and threshold behavior of the positive-energy phase shift as a 
function of the strength of the scalar square- well potential (k = — 1). 

FIG. 7. Upper (g) and lower (/) components of the lowest energy k=—1 bound state for three 
different values of the strength of the scalar square-well potential. 
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TABLES 



TABLE I. Illustration of Levinson's theorem for k = — 1 using a mean-field approximation to 



the Walecka model for 40 Ca. 
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S k (+M)/tt 




a4 +) 


6 K (-M)/ir 


a4 o) 




0.20 
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1 
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0.60 
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1.00 
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1.40 


2 
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10 


10 
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